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Abstract. We propose and study a scalar extension of the Standard Model which respects 
a Z 3 symmetry remnant of the spontaneous breaking of a global t/(l)DM symmetry. Con¬ 
sequently, this model has a natural dark matter candidate and a Goldstone boson in the 
physical spectrum. In addition, the Higgs boson properties are changed with respect to the 
Standard Model due to the mixing with a new particle. We explore regions in the parameter 
space taking into account bounds from the measured Higgs properties, dark matter direct 
detection as well as measurements of the effective number of neutrino species before recom¬ 
bination. The dark matter relic density is determined by three classes of processes: the usual 
self-annihilation, semi-annihilation and purely dark matter 3 —)• 2 processes. The latter has 
been subject of recent interest leading to the so-called ‘Strongly Interacting Massive Particle’ 
(SIMP) scenario. We show under which conditions our model can lead to a concrete realiza¬ 
tion of such scenario and study the possibility that the dark matter self-interactions could 
address the small scale structure problems. In particular, we find that in order for the SIMP 
scenario to work, the dark matter mass must be in the range 7—115 MeV, with the global 
symmetry energy breaking scale in the TeV range. 
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1 Introduction 

The existence of dark matter (DM) comprising approximately 85% of the matter content in 
the Universe, as supported by several astrophysical and cosmological observations, is a strong 
evidence that the highly successful Standard Model (SM) of particle physics is incomplete: 
there must be new physics beyond the Standard Model (BSM physics) (for a review, see 
e.g. [1-5]). 

If DM is made out of new particles, these must be electrically neutral, have a lifetime 
longer than the age of the Universe and be produced in the observed amount. There are several 
models containing particles that are suitable candidates for DM, such as supersymmetric 
(SUSY) extensions of the SM with a discrete R-parity symmetry that makes the lightest 
supersymmetric particle stable. However, at this point there is no evidence for this class of 
models from searches conducted at the Large Hadron Collider (LHC) and the simplest SUSY 
models are under stress due to naturalness issues (see, e.g. [ 6 ]). 

Another class of BSM scenarios with a DM candidate consists in extending the SM with 
an augmented scalar sector furnished with a discrete symmetry ensuring the stability of the 
DM particle. If the extra scalars are postulated to be SM singlets, they can interact with SM 
matter only through their couplings to the Higgs doublet, the so-called Higgs portal scenarios. 


- 1 - 




This class of models may lead to invisible Higgs decays, allows for strongly self-interacting 
DM [7, 8 ] and may also alleviate the vacuum stability issue of the SM [9]. 

The simplest model in this class is the SM with an extra real scalar field respecting a 
discrete Z 2 symmetry [ 10 - 12 ]. Constraints on this scenario have been recently updated [13]. 
In fact, requiring the model to respect the observed DM abundance, direct detection bounds 
and the limits on the invisible Higgs width severely constrains the parameters, placing the 
extra scalar mass in the range between 53 and 63 GeV (see also [14]). 

Another possibility is to consider scalar DM models with a Z 3 symmetry. One particular 
feature of these models is that they allow for semi-annihilation processes [15-18], where the 
annihilation of two DM particles can give rise to a DM particle in the final state, as for example 
XX AT, where X and T stand for a DM and a non-DM particle respectively. In addition, 
they predict other interesting number-changing processes involving only DM particles, such 
as XXX —)• XX. These inelastic processes can be relevant for the chemical equilibrium 
and freeze-out of DM particles and their effect must be included in the Boltzmann equation 
controlling the DM density. Remarkably, Z 3 models with that feature can also address the 
‘core versus cusp’ [19-22] and ‘too big to fail’ [23, 24] problems because they may exhibit 
sizable DM self-interactions. In this case, the DM in this type of scenarios has been dubbed 
Strongly Interacting Massive Particle (SIMP) [25] Z 3 models have been recently studied, 
assuming that this symmetry is either just a global symmetry as in [27] or that it is a remnant 
of a local 17(1)dm symmetry spontaneously broken to Z 3 [28]. 

In this paper we explore an intermediate situation, where a global 17(1)dm symmetry is 
spontaneously broken to a discrete Z 3 symmetry by postulating the existence of a scalar field 
with charge 3 under the global group and a non-vanishing vacuum expectation value [29] 
Due to this, in addition to a natural DM candidate there is a Goldstone boson (GB) in the 
physical spectrum. In our study, we include the novel processes mentioned above in the 
determination of the DM relic abundance and consider the role of GB in its production. In 
doing so we take into account bounds from the measured Higgs properties, DM direct detection 
as well as measurements of the effective number of neutrino species before recombination. In 
addition, we find under which conditions our model can be a concrete implementation of the 
SIMP scenario. 

This paper is organized as follows. In the next section we introduce the model and 
discuss the different mechanisms that could generate the DM relic density. In section 3 we 
analyze the constraints arising from the effective number of neutrino species, Higgs physics 
and direct detection experiments. Following that, in section 4 we focus on the self- and semi- 
annihilating scenarios and show that in these cases DM behaves like a Weakly Interacting 
Massive Particle (WIMP). Subsequently in section 5, we address the possibility of having DM 
as a Strongly Interacting Massive Particle (SIMP). In particular we study how to implement 
the 3 —)• 2 annihilation mechanism for its production and show that this naturally leads to self¬ 
interacting cross-sections that typically address small scale structure problems in astrophysics. 
Finally, in section 6 we present our conclusions. 

A model based on QCD-like strongly coupled gauge theories with a Wess-Zumino-Witten term was recently 
put forward in [26]. 

■^The origin of the t/(l)DM group could be attributed to an underlying flavor symmetry in order to generate 
radiative neutrino masses and to stabilize the DM as recently discussed in [30]. 
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2 Z 3 Dark Matter from the Breaking of a Global t/(l)DM 
2.1 Description of the Model 


We postulate a dark sector with a global C/(1)dm which is spontaneously broken into a Z3 
symmetry. This is achieved by considering, in addition to the usual scalar doublet H, two 
complex scalar fields (px and X which are singlets under the SM group. While the SM 
particles do not transform under the 17(1 )dM) the new scalars (px and X have charges equal 
to 3 and 1 respectively. This choice leads naturally to a Z 3 symmetry after the spontaneous 
breaking of the continuous group [29]. In fact, if we assume that in the vacuum state {X) = 0 
and {(px) = / 0 , then the breaking SU{2)l <8* U{1)y <81 ?7(1)dm — t C/(1 )em <8> Z3 occurs. 

To see this explicitly, let us consider the most general renormalizable scalar potential 

V=- p*xPx + X4P*x(Pxf + lAx*X + Xx{X*Xf 

+ X^H^x^x + X^xX*Xp\Px + XhxX*X H^H + {X^X^p*x + H.c.) . (2.1) 

After symmetry breaking we can write the fields in the unitary gauge as 



0 

VH + h{x) 


(px 


1 

V2 


[v^ + p + ir]), 


( 2 . 2 ) 


where vh = 246 GeV is the vacuum expectation value of H, which gives mass to the SM 
particles via the Brout-Englert-Higgs mechanism [31]. Therefore the part of the potential 
involving X becomes 

VA = m\X*X + Xx{X*Xf + ]^Xhx {^VHh + X*X 

X*X + 

where 

+ 2^<t>xv‘^ + -XhxVr ■ ( 2 - 4 ) 

The potential in Eq. (2.3) is manifestly invariant under the remnant Z 3 C 17(1)dm that 

27r • 

acts non-trivially only on the field X in the following way: X ^ e 3 It is precisely 
because of this reason that the particle associated to X can not decay and is thus identified 
with the DM. In other words, the Z 3 group ensures the stability of the DM in the present 
model. Notice also that while X and X belong to different representations of Z3, they have 
the same mass and therefore they both contribute to the total amount of DM. 

Due to the X^pu coupling there is a mixing between the h and p fields, which determines 
the masses for the physical states h and p 

~ v'h Xfjyv'^ + {Xh v'jj — X(j) u^)^ + {X(f,H Vh , (2-5) 

ml = XHvjj + X^v^- {Xh vjj - X^ u^)2 + (A,^^ vh , (2.6) 


+ 2^<t>x (2f</,p + 


A 3 

V2 


X^(u 0 + p - ir/) + H.c.^ , 


(2.3) 


where we identify h with the SM Higgs boson with mass ruh 
9 is defined by 


tan 29 


X<pH VH V<j} 
A<i> ~ vjj 


r\j 


125 GeV and the mixing angle 


(2.7) 
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From this it follows that the masses and are related by 


V 


2 

4 > 


m\ sin^ 0 + TTT-p cos^ 9 

2 


( 2 . 8 ) 


On the other hand, the held rj remains massless and it is hence the GB associated to the 
symmetry breaking of the global f7(l)DM- 

Consequently, the physical spectrum of this model consists of the real scalars h, p and 
the GB ? 7 , along with the complex scalar X. There are ten free parameters in the potential. 
Fixing rrih and vh to their physical values reduces the number of independent parameters to 
eight. We will use the following parameters to characterize the model: two masses {nip and 
mx), the mixing angle 9 and hve quartic couplings A^, Ax, A 3 , A^x and \hx- According 
to the potential ( 2 . 1 ), these couplings are real with the exception of A 3 , whose phase can 
be absorbed by a held redehnition, and therefore the dark sector does not introduce any 
CP-violation. 

Let us point out that, with the exception of the term involving A 3 , all the terms of 
Eq. (2.3) contain pairs of DM helds. Hence, the A 3 coupling characterizes the Z 3 DM phe¬ 
nomenology. On the other hand, Xhx controls the so-called Higgs portal because it connects 
the SM with the DM particles via the Higgs boson. Likewise, A^x connects the p and rj 
bosons with the DM, and we therefore call it the GB portal. 

Some comments are now in order. Stability of the potential requires that it should be 
bounded from below for large values of the helds. One could also require perturbativity of 
scattering amplitudes up to high energy scales, such as the Planck scale, which typically 
imply that the coupling constants should not be too large at the electroweak scale. Some of 
these issues were discussed in [28]. For the purposes of this phenomenological analysis, we 
just assume that these conditions are met; however, we do require our couplings to be in the 
perturbative regime, |A| < dvr at the electroweak scale. 

The model was implemented in FeynRules [32, 33] and the output was used both in 
CalcHEP [34] and MicrOMEGAs [35, 36] in the phenomenological studies discussed in the fol¬ 
lowing sections. 


2.2 Dark Matter Relic Abundance 

Due to the Z 3 symmetry, in addition to the usual self-annihilation of X and A, there are 
other processes that change the number of DM particles, and hence contribute to the thermal 
abundance of DM. For instance, two X particles can semi-annihilate into a X and a Z 3 singlet 
state, or three X particles can annihilate into a XX pair. All these processes can be divided 
in three classes: self-annihilation [37], semi-annihilation [15, 16] and 3 —)• 2 annihilations [25]. 
The first category corresponds to the usual processes where two DM particles annihilate into 
non-DM particles. In addition, semi-annihilation occurs when two DM particles produce a 
single DM particle in the final state. Finally, 3 —)• 2 processes involve only DM particles both 
in the initial and final states. 

If we assume no asymmetry between X and A, the total DM number density n associated 
to both particles satishes the Boltzmann equation [25, 38] 


:h+3ff„ = -2 
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13^2 , (2.9) 
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Figure 1. Freeze-out of the DM particles for the different thermal production modes, for 

self-annihilation, semi-annihilation and the 3—^2 process, giving rise to the measured relic abundance 
(see text for details). The equilibrium density is also depicted. A DM mass mx = 1 GeV was assumed. 


where H is the Hubble parameter and neq is the equilibrium number density for a given DM 
temperature T. The thermally-averaged interaction rates 7 are defined in Appendix A, and 
are given in terms of the cross-sections by 


dself — 2^6q(^^)self 1 Tsemi — ^eq(^^)semi 1 

73^2 = neq((Tu)2^3 = V^)■ 


( 2 . 10 ) 


The three terms on the right-hand-side of Eq. (2.9) correspond to the self-annihilation, 
the semi-annihilation, and the 3 —>• 2 annihilation processes, respectively. Here {av) is the 
standard thermally-averaged velocity-weighted cross-section for 2 —)• n reactions and {(TV^) 
is the equivalent quantity for 3 —)• n processes, as discussed in Appendix A. Let us note 
that the Boltzmann equation should contain other terms corresponding to three DM particles 
annihilating into a two-body final state with none or one DM particle. However these inter¬ 
actions are subdominant (c.f. Appendix B) and will be neglected in the following analysis. 
Notice that {av) 2^3 = neq {av‘^) 3^2 as shown in Appendix A using CP conservation and the 
principle of the detailed balance [39]. 

In this work we solve the Boltzmann equation (2.9) in order to estimate the DM relic 
density abundance today, which is given by 

Hdm = (2.742 • 10® GeV^) - my , (2.11) 

s 

where s is the entropy density of the photon plasma. That value must agree with the one 
measured by the Planck collaboration [40]: 

= 0.1196 ±0.0031. (2.12) 


For the sake of illustration, Fig. 1 shows the evolution of the comoving number densities 
of the DM as a function of x = mx/T obtained by solving the Boltzmann equation for 
the different thermal production modes considered separately: self-annihilation (red), semi¬ 
annihilation (blue) and the 3 —)• 2 process (orange). We take mx = 1 GeV with constant 
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Figure 2. Left Panel. GB contribution to the effective number of neutrino species as a function 
of its decoupling temperature from the thermal plasma. The dashed vertical line corresponds to 
= m^. Right Panel. Mixing angle as a function of the p boson mass. The continuous lines 
represent T^ = The dashed lines correspond to the upper limit given by the invisible decay of 
the Higgs, assuming mx > rn-hjl so that the Higgs and the p can not decay into a couple of DM 
particles. In addition, the hatched region shows the excluded mixing angles from p searches in the 
OPAL detector at LEP. 

{av)seii, (o'f)semi and chosen so that the final abundance equals the measured one. 

Since this plot is only for illustration, we assume that all the relativistic particles have the same 
temperature (this might not always be the case, as shown in section 3.1). The corresponding 
equilibrium density is also depicted in black. As shown in the figure, although the 3 —>• 2 
process freezes out at larger x compared to the self- and semi-annihilations, the former reaches 
its final DM abundance much faster. 

In sections 4 and 5, we analyze further the different regions of parameter space where 
the self-annihilation, semi-annihilation and 3 —)• 2 processes are dominant and study the 
phenomenological consequences, taking into account the constraints discussed in the following 
section. 

3 Constraints on the Model 

In this section we discuss some constraints arising from the production and decay of the SM- 
like Higgs, DM direct detection, as well as from the measurement of the effective number of 
neutrinos species before recombination. 

3.1 Effective Number of Neutrino Species 

Being massless particles, GBs can fake the effects of neutrinos in the Cosmic Microwave 
Background (CMB). The existence of these new relativistic degrees of freedom, the so-called 
Dark Radiation, is characterized by the effective number of neutrino species N^s (for a recent 
discussion see, e.g. [41]). Measurements of the anisotropies in the CMB performed by the 
Planck collaboration in 2013 indicate that N^s = 3.30 ± 0.27 [40], consistent with the SM 
prediction In this model, the contribution to the value of A^eff from the GB is determined 
by the moment in which it goes out of equilibrium with the thermal plasma If and Tij 

^After the completion of this work, new preliminary results form the Planck collaboration indicate an even 
better agreement with the SM prediction A7ff = 3.13 ± 0.32 [42]. These results are consistent with BBN [43]. 

^The equilibrium between the SM and the dark sector is achieved in the Early Universe by a number of 
processes. For instance by hh —>■ hp, which is only suppressed by the mixing angle 6, and which could establish 
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Figure 3. Diagrams for the GB annihilation and scattering relevant for the thermal decoupling. 


are respectively the decoupling temperatures of the GBs and neutrinos from the plasma and 
if the GBs are not reheated after their decoupling from the SM, then 


AiVeff 


Aeff - 3 = 


4 

7 


'giSl 

.9iTp^) 



(3.1) 


where g{T) is the number of relativistic degrees of freedom at a given temperature. We plot 
the resulting contribution as function of the decoupling temperature on the left panel of Fig. 2. 
From this it is clear that when the decoupling temperature is greater than the muon mass, 
A^eff is reasonably close to the SM prediction. Gonsequently, we have to make sure that the 
interaction rates keeping the equilibrium between the GB and the SM become smaller than 
the Hubble expansion rate at temperatures greater than the muon mass. In fact, assuming 
nip greater than a few GeV is enough to have a diluted GB contribution to N^s. 

In order to show this, we consider the case ~ nip (For a detailed analysis, see e.g., 
[44, 45]). Then the equilibrium takes place via reactions such as r] f rj f and // 0 gr], 
where / denotes a SM fermion, as shown in Fig. 3. Since the corresponding amplitudes are 
proportional to the fermion mass, processes involving muons are dominant. Then, the depar¬ 
ture of equilibrium occurs when the Hubble expansion becomes comparable to the interaction 
rate with muons. Exploiting the fact that the latter depends only on the three parameters 
nip, Xff, and 9 (as well as on T^), in Fig. 2 we show the values of the mixing angle for which 
the interaction rate equals the Hubble parameter for different values of A^. The regions below 
the solid lines correspond to decoupling temperatures higher than the mass of the muon and 
thus to a diluted contribution to A^eff- From Fig. 2 we conclude that such contribution is 
small when the p scalar is heavier than about 5 GeV or when the mixing angle is smaller than 
about 10“®. 


3.2 Higgs Sector 

The enlarged scalar sector increases the number of Higgs decay channels. The new decay 
modes are into pairs of p, g and DM particles, when kinematically allowed. The last two 
channels [h ^ gg and h ^ X X) contribute to the invisible decay of the Higgs, whereas the 
first one (/i —)• p p) only does it partially when the p pair itself decays into a couple of g or 

the equilibrium at temperatures above the TeV scale. We corroborate this for the scenarios considered in this 
work. 
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Figure 4. Tree level diagrams for the decay h XXX. 


DM particles. The corresponding decay rates are: 


T{h ^ = 


327r 


r(/i -)■ pp) = — 2 2^ 2 cosO - sin^)^ sin^ 26 , 

i.Zo7r ^ ^ 


T{h 


J m\ — Am\ 2 

^ on _,...,2 - Vcj) sin 6 - Xhx vh cos Of 


327r 


(3.2) 

(3.3) 

(3.4) 


Notice that three DM particles form a Zs-singlet and as a consequence the decays h —>• XXX 
and h —>• XXX are allowed. The corresponding Feynman diagrams are shown in Fig. 4, 
where one can see that its rate is proportional to A 3 . Nevertheless, this channel is phase- 
space suppressed and therefore subdominant compared to the two-body decay modes. 

Correspondingly, the decay channels for the p boson are 


r(p^ 7777) 

T{p^ XX) 


mp cos^ 6 

327r 77? 


mp — 4m‘^ 

32tt mp 


{^<px cos 6 -h Xhx vh sin 6)" 


(3.5) 

(3.6) 


In addition, due to the h — p mixing, the p boson can also decay into SM particles, 
with a rate give by r(p —>• SM SM) = T{hsM SM SM) x sin^ 6, taking mhgy^ —>• mp in the 
right-hand-side of the expression. Because of the same reason, all the decay widths of the h 
boson into SM particles decrease compared to the SM ones by a factor cos^ 6. 

With these decay rates, we calculate the invisible Higgs branching ratio and apply the 
bound reported in [46] . That analysis considers simultaneously a universal modihcation of the 
Higgs boson couplings to SM particles as well as the possibility of an invisible Higgs decay, 
concluding that 

[1 - BRinv] > 0.81 (@ 95 % C.L.), ( 3 . 7 ) 

where k is the universal modification for the Higgs coupling. In our case we identify k = cos 9 
and the invisible Higgs decay branching ratio with 

BRinv = BR{h ^ 77 77) BR(/i ^ X X) + BR(/7 ^ X X X) + BR{h ^ X X X) 

+BR(/7 ^ pp)x [BR(p ^ 7777) + BR(p ^XX)]^ . 
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Figure 5. Diagrams responsible for DM direct detection. 


Although this is the most stringent constraint, we also require that the total decay width 
never exceeds the experimental bound [47] 

< 22 MeV (@ 95% C.L.). (3.9) 

Direct searches for a light neutral scalar produced in association with a Z in the OPAL 
detector [48] at LEP also apply. By identifying such a particle with the p scalar, it is possible 
to set an upper bound on the mixing angle. 

We apply these limits for all the scenarios considered in this work. In the right panel 
of Fig. 2 we show them in the particular case when the scalars can not decay into DM 
particles. The dashed lines represent the upper bounds arising from the Higgs invisible decay 
of Eq. (3.7), whereas the upper hatched region is excluded because of direct searches of the p 
scalar in OPAL. 


3.3 Dark Matter Direct Detection 


The scattering of DM particles off nuclei in direct detection experiments takes place via the 
t-channel exchange of h and p bosons as depicted in Fig. 5. Such process only generates a 
spin-independent signal with a DM-nucleon cross-section given by 


O'XN 


f%m% + m?p) + {\hx u//cos 20 - \^xv<psin2e){m?p - m|)]' 

dvr v'jj {mx + rnx)^ 


(3.10) 


where mx denotes the nucleon mass and fx ~ 0.27 is a constant that depends on the 
nucleon matrix element [36]. For GeV Dark Matter, the resulting expression is constrained 
by comparing it to the upper bound obtained by the LUX collaboration [49]. 


4 Weakly Interacting Dark Matter 

If self- or semi-annihilations dominate the production of DM, the latter is in kinetic equili¬ 
brium with the SM particles and the GB (See Appendix A). As discussed in section 3.1, due 
to the constraint on this equilibrium can only exist at temperatures greater than the 
muon mass. Since the freeze-out temperature is roughly mx/25, we thus consider only self- 
and semi-annihilations in the GeV range. 

Because of this mass range and since the couplings involved in DM phenomenology 
are typically small, self and semi-annihilation scenarios describe Weakly Interacting Massive 
Particles (WIMP). We now discuss them separately. 

4.1 Self-annihilating Scenario 

If the self-annihilation of DM dominates over the other two processes, that is if 7seif is much 
greater than 7semi and 73 _). 2 , the DM production proceeds via the familiar Lee-Weinberg sce¬ 
nario [37]. In this case, as shown in Fig. 6, the DM abundance is obtained by self-annihilation 
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Figure 6. Tree level diagrams for the annihilation between X and X. Here f corresponds 
to any SM fermion and V to any massive gauge boson. 


into pp, pp, hp and pairs of SM particles, until the reactions freeze-out at x ~ 25. DM an¬ 
nihilates into SM fermions and vector bosons via the s-channel exchange of a /i or a p boson. 
For the other annihilation channels the contact term and the t-channel exchange of a X are 
also present. The only exception is the reaction XX ^ pp where the t-channel exchange 
does not exist. All these diagrams are proportional to Xhx and X^x as opposed to the ones 
corresponding to the semi-annihilation and the 3 —)• 2 mechanisms which are proportional to 
A 3 . As a result, self-annihilation dominates whenever A 3 is much smaller than Xhx and X^x- 

If we neglect the mixing angle, there are two regimes depending on the relative size 
of the Higgs and GB portals Xhx and X(j)X- On the one hand, when the latter dominates, 
DM annihilates into the dark sector, that is, into GB and p pairs. The two channels are 
democratically distributed provided that mx S> mp. On the other hand, when Xhx S> X^Xj 
the annihilation produces SM particles pairs. Due to the s-channel Higgs exchange, the 
heaviest kinematically allowed pair tend to be produced. These two regimes are illustrated in 
the upper panel of Fig. 7, where we show the parameter space giving rise to the measured DM 
relic abundance, for different combinations of mx, Xhx and A^x- There we take mp = 5 GeV, 
sin0 = 10 “^, X^ = 0.1 as well as A 3 = 0 in order to avoid semi-annihilations and 3 —)■ 2 
annihilations. We calculate the relic abundance using MicrOMEGAs-3. From the figure it is 
possible to see that large masses in general require large couplings. We can also see that when 
Xhx dominates, due to the existence of the Higgs resonance and kinematic thresholds, the 
dependence on the mass is more complicated than in the other case. This is exemplified in the 
lower panels of Fig. 7, where the dependence on the DM mass is shown explicitly. The Higgs 
boson funnel appears as well as the kinematic openings of the annihilations into W~^ W ~, 
Z Z and hh. The LUX experiment rules out parts of the parameter space corresponding to 
DM masses O ( 10 ) GeV and sizable Xhx couplings. This corresponds to the hatched regions 
in Fig. 7. 

Similarly, the left (right) panel of Fig. 8 shows in blue the regions of the parameter space 
that could give rise to the measures DM relic abundance in the plane \mx, Xhx] {[np^x, 
after marginalizing with respect to X^x (Xhx)', that is, each plot includes all the possible 
values of the corresponding quartic coupling. The empty regions, corresponding to high values 
for the couplings, always generate a DM under abundance. In addition, both the LUX and the 
Higgs invisible decay exclusions are shown in the figure as hatched regions. Again, those areas 
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Figure 7. Self-annihilation scenario. Parameter space giving rise to the measured DM relic 
abundance, for different combinations of mx and the Higgs and GB portals Xhx ^nd A^x ■ We have 
chosen mp = 5 GeV, sin0 = 10“^ and = 0.1, we also set A 3 = 0 in order to avoid semi-annihilations 
and 3 —>■ 2 reactions. The hatched regions are ruled out by the LUX experiment. 


correspond to the points that could be excluded for some combinations of the parameters. 
The unhatched regions are always allowed. 

4.2 Semi-annihilating Scenario 

In the previous section we assumed that A 3 was negligible with respect to Xhx and X^px- 
In general this may not be the case and therefore the semi-annihilation and the 3 —>• 2 
process have to be taken into account. Although the latter process can still dominate, here 
we simply assume that it is subdominant and work under that hypothesis. In section 5.2 
we will address the general case and discuss under which circumstances the 3 —>• 2 reaction 
can be safely ignored. Under these conditions, the generation of relic abundance occurs not 
only via the self-annihilation of DM particles but also via the semi-annihilation processes 
DM DM —)• DM S, where S could be h, p or p. The corresponding Feynman diagrams at 
tree-level are shown in Fig. 9. Notice that the semi-annihilation into GBs is always possible 
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Figure 8 . Self-annihilation scenario. Same as Fig. 7 hut marginalizing with respect to X^x 
(^HX) in the left (right) panel. We also include the constraint from the Higgs invisible decay width. 
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Figure 9. Tree level diagrams for the semi-annihilation process X X ^ X S where S could 
be h, p or ij. 


and can only be suppressed by taking A 3 very small. 

When A 3 is much greater than X^x and Xpx, self-annihilations can be neglected as 
well as the second and third diagrams of the semi-annihilation process in Fig. 9; in other 
words, only the contact interaction diagram contributes. Hence, the relic abundance scales 
like ~ l/((Tu)semi ~ {^x/^ 3 )'^■ This behavior is illustrated in Fig. 10, where we depict 

the regions of the parameter space [mx, A 3 ] giving rise to the observed relic density when 
X(f,x = ^HX = 0. There we also show the equivalent regions when we relax this assumption, 
while still keeping A^x = Xhx for the sake of simplicity. In this plot we choose = 5 GeV, 
sin0 = 10“^ and X^ = 0.1. In the case where X^x = ^HX are sizable, both the Higgs boson 
funnel and the kinematic openings of the annihilations into W~^W~, ZZ and hh are visible. 
Similarly, the LUX experiment rules out parts of the parameter space corresponding to DM 
masses 0(10) GeV. These parts are shown in the Fig. 10 as hatched regions. 

As we did for the self-annihilation case, in Fig. 11 we show the region of the parameter 
space that could give rise to the measures DM relic abundance in the plane [mx, A 3 ] after 
marginalizing with respect to Xpx and A(^x- The white region, corresponding to large values 
for A 3 , always generates a DM underabundance. Moreover, assuming that the cross-section 
for the semi-annihilation into a GB is never greater than about (cru)semi ~ 3 • 10“^® cm^/s, 
this translates into the upper bound 


A 3 < 3 • 10-^ 


mx \ 
IGeV/ 


(4.1) 


which can be seen in Fig. 11. In addition, we also show both the LUX and the Higgs invisible 


12 


























Figure 10. Semi-annihilation scenario. Parameter space giving rise to the measured DM relic 
abundance in the plane [mx, A 3 ] for different combinations of the couplings = ^hx- We have 
chosen mp = 5 GeV, sind = 10“^ and = 0.1. The hatched regions are ruled out by the LUX 
experiment. 



Figure 11. Semi-annihilation scenario. Same as Fig. 10 but marginalizing with respect to X^x 
and Xhx ■ We also include the constraint from the Higgs invisible decay width. 


decay exclusion limits. 

5 Strongly Interacting Dark Matter 
5.1 Dark Matter Self-interactions 

In spite of the great success of the ACDM scenario in explaining the formation and evolution 
of cosmic structures, there are some problems at small scales (see, e.g. [50]). Simulations of 
this scenario result in cuspy density profile of halos towards its center and a large number of 
small satellite halos. These results are challenged by observations. Although the inclusion of 
baryonic effects can ameliorate the discrepancies, it is still unclear whether it is necessary to 
change the physics of the DM sector. For instance, self-interacting DM with a strength [Sl¬ 
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Figure 12. Diagrams for the elastic scattering process X X ^ X X and XX X X. 


58]: 

0.1 cmVg < < 10 cmVg (5.1) 

obs 

can help in solving the observed discrepancies. Yet, the Bullet Cluster [59-61] and recent 
analysis of the constraints from halo shapes [54-56], favor the following upper bound on the 
DM self-interacting cross section (at velocities greater than 300 km/s) 


^scatter 

mx 


< 1 cmVg. 


(5.2) 


The corresponding elastic scattering processes in our model are shown in Fig. 12. From 
the diagrams it is possible to see that the exchange of a p or a h boson lighter than the DM 
can possibly lead to non-perturbative effects such as long range interactions. In this section 
we limit ourselves to the perturbative scenario and we therefore assume mp, mu mx- A 
posteriori we will find that this is also necessary for the DM production to take place via the 
3 —>• 2 annihilation process (see section 5.2). Under this assumption the elastic scattering 
cross-section in the non-relativistic limit reads 


•^scatter — 


3 (4 -20 Ax + 

647r 


(5.3) 


where Z = \^v^/mx- The Z parameter can not be arbitrarily large. In fact, perturbative 
unitarity requires the s-wave amplitude of the process XX —)• XX to be bounded from 
above [62]. Considering rup and mx, this translate into 


’ A/vx VH sin e + X^x cos 9 ' 

2 

’ A/vx VH cos 9 - X^x V^ sm9' 


“T 

mh 


2Xx + 9Z^ + 


which implies the following weaker bounds 


< dvr, (5.4) 


Z < — ~ 1.2 and Ax < 27r . 

3 


(5.5) 


The process XX —>• XX also sets an upper bound on the same parameters, however it is 
never as stringent as the one of Eq. (5.4). 
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Figure 13. Self-interacting cross-section. Regions of the parameter space preferred by the 
astrophysical observations of small scale structures anomalies, for different values of the DM mass 
(see Eq. (5.1)j. The areas disfavored by the Bullet Cluster (i.e. Eq. (5.2) j are also shown. In the 
hatched region perturbative unitarity is lost. 


Fig. 13 shows the regions of the plane [Ax, -2^] fulfilling Eq. (5.1) for different values of 
the DM mass. In addition, we depict the areas disfavored by the Bullet Cluster for those 
masses (i.e. Eq. (5.2)) and we also hatch the regions where perturbative unitarity is lost. In 
this figure a limit on the DM mass is manifest, corresponding to 

mx<173MeV, (5.6) 

which is the mass that saturates the perturbative unitary bound of (5.5) and that is still in 
agreement with the astrophysical lower bound of Eq. (5.1). 


5.2 3 —>• 2 Dark Matter Scenario 


Instead of the usual 2 —)■ 2 annihilation scenarios described in previous sections, another 
mechanism for the reduction of the comoving DM number density is the annihilation via 
3 —)• 2 processes. Here we discuss this case, assuming that it is the dominant process for 
the relic density computation, or equivalently that 7seif, 7semi ^ 73 -i .2 at the epoch of the 
freeze-out. 

As explained in the Appendix A, it is necessary for the DM to release the kinetic energy 
produced in the 3 —>• 2 reactions. This can be fulfilled if DM is in kinetic equilibrium with 
either SM particles or the GBs. The first possibility is realized by means of DM scattering 
off electrons. However, it requires a large Higgs portal. This can be understood from the 
effective Lagrangian between electrons and DM. Assuming mh, rup S> mx and 0 <C 1, we 
have 


C 


^HX rne 


m 


2 

h 


X*Xee, 


(5.7) 


and one can write 


(av] 


kin — 


with e = 


\hx me mx 


m 


X 


y/^m 


(5.8) 


where me is the electron mass. In order to keep the kinetic equilibrium between the DM and 
the electrons from this interaction, it is necessary that e > 5 x 10~® [25], or equivalently, that 
Xhx ^ 1-5 X 10^ (me(mx), which is too large if Eq. (5.6) holds. 
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Figure 14. Tree level diagrams for the annihilation XX X X X . The first two rows 
contain processes that are mediated via the exchange of h, p and X whereas the last two rows involve 
solely DM particles. 


As a result, the only possibility is to have kinetic equilibrium between the DM with the 
GBs. The corresponding rate depends on X^px and for m^, mp 3> mx is mostly driven by 
a contact interaction (see Fig 6). Because of this, it is not suppressed by a heavy scale in 
contrast to the previous case. Moreover, although this scattering rate should be greater than 
the Hubble parameter at the epoch of the freeze-out, we must insure that the rate for the 
self-annihilation process into GBs -which is obtained by crossing symmetry- is negligible. This 
is possible because rkin/Tann ~ np/nytM 1 at the freeze-out, due to the sub-dominance of 
riDM during the radiation-dominated epoch. We conclude that Xpjx must be sizable but not 
too large in order to avoid self-annihilations into GBs. 

In section 3.1, it was implicitly assumed that the GB did not heat up after its decoupling 
from the SM plasma. In fact, this is not the case here because DM is not in kinetic equilibrium 
with the SM and transfers all its entropy exclusively to the GBs when it annihilates via the 
3 —)• 2 mechanism. As a result, Eq. (3.1) must be modified accordingly in order to account 
for this effect and thus T!^ must be larger than the muon mass so that the contribution to 
Aeff is negligeable. This always takes place here because we consider mp S> mx ~ 
100 MeV, which implies that the processes establishing the equilibrium are very suppressed 
for temperatures around the muon mass. 

With this in mind, we can consider the 3 —)• 2 reactions in this model. All of them 
are equivalent, up to CP conjugation, to either XXX ^ XX or XXX ^ XX. In 
Figs. 14 and 15 we show the corresponding diagrams at tree-level. The corresponding cross- 
sections can be calculated analytically if we work in the non-relativistic approximation (see 
Appendix B). In fact, the resulting expressions can be further simplified in the case 7seif ^ 
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Figure 15. Tree level diagrams for the annihilation XX X X X. 


73_>.2 because of the following reason. As for self-annihilation, some of the 3 —)• 2 diagrams 
are proportional to Xhx and X^x- In fact, this happens for those in the first and second rows 
of Figs. 14 and 15 simply because they involve exchange of p or /i particles. In contrast, the 
diagrams of the third and fourth lines depend uniquely on different couplings, namely A 3 and 
Ax- Consequently, the contribution of the first set of diagrams is presumably negligible, if 
self-annihilations are subdominant. The fact that and S> mx reinforces this argument. 
However, notice that for this to be true it is not necessary that A(^x vanishes. Under these 
assumptions we find 


{av^ 


' 3^2 — 


TbSvr 


4 (2Ax + 9 ^ 2 )" + ^ (74Ax - mZ^) 


2\2 


+ O {Xhx , Xrpx)'^ , (5.9) 


where Z = X^v^/mx was already introduced. Because self- and semi-annihilations into 
GBs are always kinematically allowed, Eq. (5.9) must be compared with the cross-sections 
associated to these processes 


(o'u)seif = O {Xhx , X^xf' , 


((Tu)semi 


27 Aj 
2567 r m\ 


(5.10) 


We already discussed the first process. Regarding the second one, although (<Tr;)semi and 
(cJU^) 3->.2 have different mass dimensions and can not be compared directly, it is clear that 
the semi-annihilations are suppressed for small A 3 , and that the 3 —>• 2 is important when Z 
is large. Moreover, a suppression in the 3 —?> 2 process due to A 3 can be compensated with a 
large v^f) with respect to the DM mass. Under these circumstances we can make sure that 7 seif, 
7 semi ^ 73 ^ 2 ) while still having kinetic equilibrium between the GBs and the DM. We remark 
that the DM mass is not required to be in the GeV range because we assume no equilibrium 
between the GBs and SM plasma and consequently there is no significant contribution to Aeff- 
In order to study the regions of the parameter space where the 3 —)• 2 reproduces the 
observed relic abundance, instead of solving numerically the Boltzmann equations directly we 
work in the freeze-out approximation. 
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5.3 The 3 —2 Preeze-out Approximation 

When the temperature of the plasma is much greater than the DM mass, the Boltzmann 
Eq. (2.9) has the solution n = TT-eq, which corresponds to the chemical equilibrium. However, 
once the temperature drops below the DM mass, the rate for the process 2 —?> 3 is kinematically 
suppressed compared to the 3 —)• 2 reaction. As a result the DM number density deviates from 
its equilibrium value. In fact, it eventually becomes much greater than its equilibrium value 
n ;:?> rieq. Hence, for these temperatures and when there are no self- and semi-annihilations, 
Eq. (2.9) can be approximated by 

^ + 3Hn^ -n^{av‘^)3^2- (5.11) 

When the DM particles are non-relativistic, the cross-section (cru^) is independent of the 
temperature. In that case, using standard methods, Eq. (5.11) admits the following solution 

{av'^)3^2 = (8.65GeV"®) (iS^) ’ 

where FO stands for the freeze-out, g are the number of relativistic degrees of freedom and 
X = mxiT. Here T is the DM temperature which equals that of the GB plasma. It is 
remarkable that, in contrast to the standard self-annihilation scenario, here the cross-section 
that matches the observed relic density depends on the DM mass. In order to estimate xpo; 
it is necessary to establish when the annihilation rate per particle n^^{av) 3^2 drops below the 
expansion rate of the Universe. Using Eq. (5.12) it is found that this happens when freeze-out 
temperature satisfies 


xpo = 20.6 -b log 


mx 


100 MeV 


gpo / xpQ X 1-5 

10.75/ V20.6/ 


(5.13) 


Because the second term of this equation depends logarithmically on the DM mass, it is a 
very good approximation to take xpo ~ 20 . 

It is somewhat interesting to calculate the value of {av) 2^3 associated to Eq. (5.12). To 
that end, notice that at the freeze-out the expansion rate is H = Ueq {av) 2^3 = nl^{av^) 3^2 
and as a result neq = H/{av‘^) 3 ^ 2 - According to Eq. (A.9) we obtain 

{av) 2^3 = \/H {av‘^) 3^2 = (1-08 x 10"® GeV"^) xpo 9 fo^ 

Ri 5.9 X 10"^® cm^/s. (5.14) 


In contrast to Eq. (5.12), notice that this equation is independent of the DM mass and is 
similar to the usual condition for the WIMP scenario. 


We use this formalism in order to calculate the relic density. In Fig. 16 we show the 
regions of the plane [Ax, 2] in agreement with the observed relic density for different values 
of the DM mass when the 3 —)■ 2 process dominates the DM production. They are obtained 
using Eqs. (5.9) and (5.12), taking into account an uncertainty of 25% in order to account 
for the error in the freeze-out and the non-relativistic approximations. As apparent from the 
figure, there are two regimes. On the one hand, for small values of Ax the DM abundance is 
set by Z. On the other hand, when Ax is of order one both variables are relevant for the DM 
dynamics. In particular notice that Z can not vanish. Fig. 16 also includes the region where 
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Figure 16. Freeze-out approximation for the 3 —>■ 2 scenario. Parameter space giving rise 
to the measured DM relic abundance for different DM masses. In the hatched region perturbative 
unitarity is lost. 


perturbative unitarity is lost, which implies an upper limit on the DM mass as shown in the 
figure: 

mx<115MeV. (5.15) 

Let us note that this limit, obtained from the relic density, is comparable to the one in 
Eq. (5.6) coming from demanding significant DM self-interactions. 

5.4 Self-interactions and the 3^2 Mechanism 

It is interesting that the variables that control DM self-interactions also determine the 3 —)• 2 
process, as shown in Eqs. (5.3) and (5.9). This suggests that both phenomena have a common 
origin, which we attribute to the fact that DM is strongly interacting. Fig. 17 illustrates this 
by showing the regions of the parameter space where the relic abundance is generated via 
the 3 —)• 2 mechanism and where the self-interaction cross-section simultaneously fulfills the 
astrophysical hints from small scale structures, as in Eq. (5.1). In addition, we show the 
areas disfavored by the Bullet Cluster {i.e. Eq. (5.2)). All this has been done under the 
assumption that perturbative unitarity is valid. The figure conclusively shows that the DM 
mass is bounded both from above and below 

7 MeV < mx < 115 MeV, (5.16) 

which is in agreement with the findings of [25]. Similarly, for the Z parameter we find 

O.35<Z<^0F. (5.17) 

We illustrate this scenario with three benchmark points which are defined in Table 1 and 
shown in Fig. 17. They correspond to different masses in the MeV range, which reproduce 
simultaneously the observed relic abundance via the 3 —>• 2 mechanism and the hints on the 
self-interaction cross-section. Whereas points B and C give rise to a crgcatter/^rix in agreement 
with the Bullet Cluster observations, point A is disfavored due to its small mass. In every 
case As, Xpx and X^px are chosen small enough in order to ensure that semi-annihilations 
and self-annihilations are subdominant (see Table 1). Nevertheless, we require that DM is in 
kinetic equilibrium with the GBs via the X^x coupling, as argued before. We again find that 
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Figure 17. Self-interacting dark matter and the 3 —)• 2 mechanism. Parameter space giving 
rise to the measured DM relic abundance, preferred by the astrophysical observations of small scale 
structures anomalies and satisfying perturbative unitarity. The three benchmark points defined in 
Table 1 are also shown. 

the scale at which the [/(1)dm symmetry breaks down has to be at least in the TeV range so 
that is of order one. We would like to remark that all these points are in agreement with 
collider constraints because the p boson is very heavy and the invisible decay width of the 
Higgs is very small. In fact, the most important contribution to the latter comes from the 
process h —^ XX with a branching ratio of approximately 1.5 • 10“®. The other channels {pp, 
rjr], XXX and XXX) contribute less than ~ 10“^^. 

6 Conclusions 

We have presented a scalar extension of the Standard Model (SM) consisting of two additional 
complex fields with charges one and three under a global 17(1)dm symmetry. The nonzero 
vacuum expectation value of the latter held spontaneously breaks the global symmetry down 
to a remnant Z 3 . This leads to the appearance of a Goldstone boson (GB), two Higgs-like 
particles, one of which is identihed with the SM scalar, and a complex held that, due to the 
discrete symmetry, can not decay and is therefore a candidate for dark matter (DM). 

This model is constrained by different observations. In hrst place, the enlargement of 
the scalar sector could potentially modify the dynamics of the recently found SM scalar. In 
particular, its decay into invisible particles due to opening of new channels involving the GB 
and the DM. Secondly, the exchange of the Higgs-like particles can mediate elastic scattering 
between the DM and SM particles leading to detectable signals in direct detection experiments 
for DM. Lastly, GB can mimic neutrinos in the Cosmic Microwave Background, signihcantly 
altering the so-called effective number of neutrinos A^eff- We have found that these constraints 
are always satished if the Higgs-like particles are heavier than ~ 5 GeV, and if their mixing 
angle is sufficiently small, namely smaller than about 10 “®. 

In these regions of the parameter space we have investigated the thermal production 
of DM. On the one hand, DM can self-annihilate either into SM particles by means of the 
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A 

B 

C 

mx [MeV] 

10 

50 

100 

-^3 

5 • 10 -'^ 

2 • 10 -® 

5 • 10 -® 


0.8 

10-2 

10-3 

^(p 

1 

o 

10 -^ 

10 -^ 

Z 

0.5 

0.9 

1.1 

rrip [GeV] 

141 

318 

311 

[TeV] 

10 

22.5 

22 

{(tv) 2^3 [cm^/s] 

6.6 • 10-26 

(M 

1 

o 

r—1 

00 

4.5 • 10-26 

Tsemi/ 'TS — >2 

0.03 

0.03 

0.05 

7self/ 73-)-2 

CO 

1 

O 

r—1 

1 

O 

r—1 

CN 

4 • 10-3 

O-scatter /mx [cm^/g] 

6.9 

0.91 

0.27 


Table 1. Benchmark points that generate the measured DM relic abundance mainly via the 
3—^2 processes. Here sin0 = 10“^ and Xhx = 10“"^. Similarly, in order to keep the DM in kinetic 
equilibrium until its freeze-out we set = 10“®. mp and are derived quantities. 


Higgs portal, or into the non-DM additional scalars via the GB portal. Similarly, DM can 
semi-annihilate, that is two DM particles can be converted into another DM particle and a 
non-DM scalar. Because the GB is massless, both annihilation processes into GBs are always 
kinematically open and they therefore play an important role in the DM production. For DM 
in the GeV range, the typical values for the Higgs and GB portals as well as the coupling 
responsible for semi-annihilations are similar to the weak couplings of the SM. Because of 
this reason, when one of these two processes dominates, DM behaves as a weakly interacting 
massive particle (WIMP). 

On the other hand, our model provides for a novel DM production mechanism, the 
so-called 3 —)• 2 mechanism consisting of the annihilation of three DM particles into two of 
them, which is possible due to the Z 3 symmetry. We find regions in the parameter where 
this mechanism is the dominant one. By studying the freeze-out approximation of the corre¬ 
sponding Boltzmann equation, and assuming the absence of non-perturbative effects, we have 
found that only DM masses in the MeV range can reproduce the observed relic abundance. 
In addition, for this mechanism to dominate over the other production modes, the Higgs 
and GB portals are required to be suppressed. Also, it is found that the C/(1)dm symmetry 
must break down at least at the TeV scale. This naturally leads to relatively strong cubic 
DM self-interactions. Because of this, this scenario is a realization of the strongly interact¬ 
ing massive particle (SIMP) paradigm. Remarkably, the same mass range and similar DM 
cubic couplings naturally give rise to self-interacting cross-sections which could solve the dis¬ 
crepancies between observations and simulations of small-scale structures in the Universe, as 
shown in Eq. (5.1). In fact, by performing a dedicated analysis of the parameter space we 
have found that in this scenario the DM mass must be in between 7 and 115 MeV. In order 
to illustrate our findings, we have studied three benchmark points (see Table 1), compatible 
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with all experimental constraints, in which the relic abundance is produced via the 3 —)• 2 
mechanism and where the DM self-interactions are in the range suggested by astrophysics. 

Direct and indirect detection experiments as well as collider searches will continue closing 
in on the parameter space of the WIMP scenario, due to the comparatively large couplings 
with the SM particles. Although this is more challenging for the SIMP scenario, astrophysics 
opens up a new window on the dynamics of DM self-interactions. This is the case, for instance, 
of the benchmark point A which is in tension with observations of the Bullet Cluster and 
recent analysis of the constraints of halo shapes. Consequently, new astrophysical observations 
can shed light on the SIMP scenario, potentially constraining further its parameter space. 
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A The Boltzmann Equation 


In this Appendix we consider the Boltzmann equation for the DM particle and how to calculate 
the corresponding annihilation rates. Since we consider no asymmetry between X and X, we 
can think of the Z 3 charge as an internal degree of freedom for the DM. As a result the DM 
density is given by 


n{t) = 2 J (A.l) 

where f{E,t) is the DM phase-space distribution. The evolution of the DM density is given 
by the Boltzmann equation 

+ ‘iHn = (A. 2 ) 


where H is the Hubble parameter, while and Np are the the number of DM particles in 
the initial and hnal states. In addition, the interaction rates are given by 

^■>-0 = / (n (n 24 (w) - w) X \M^X‘ . (A-3) 


where is the square invariant amplitude including possible symmetry factors in both 

the initial and final state, respectively. Without loss of generality we assume > 1. For 
instance, in self-annihilation processes = 2 and = 0, in semi-annihilations A„ = 2 and 
Nj^ = 1, and for the 3 —)• 2 reactions Na = 3 and Np = 2. 

We assume now that the non-DM particles are in chemical equilibrium at a common 
temperature T during the period of interest. Their phase-space distribution is thus given by 
yeq(^, ^ exp(—A/T). Although this is not true for the DM particles, they are assumed 
to be in kinetic equilibrium, as discussed below. If this is the case, because of symmetry con¬ 
siderations in a Friedmann-Lemaitre-Robertson-Walker Universe, their phase-space density 
f{E, t) is proportional to the chemical equilibrium distribution f^^{E, T), with a proportion¬ 
ality factor independent of the three-momentum [63]. This assumption implies 





Nc 

7q!—) 


(A.4) 


where corresponds to the rate of Eq. (A.3) when /j = Furthermore, due to the 

principle of detailed balance and CP conservation = 'yi3^a- Then, after adding over 

all the possible self-, semi- and 3 —)• 2 annihilation processes, the Boltzmann equation can be 
cast as 


^ + 3Hn = -2 

at 


r / \ 2 1 



/ \ 2 



n 


f n \ 

n 

— - 1 

7self - 

n ^ 

7semi — - 

-1 

\neqj 

^eq 

T-eq 

\neqj 

X'eq 


73^2 (A.5) 


or as 


dn 


— + 3Hn = -((Tu)self (n^ - ?T-eq) “ {crv) serai {n^ - nHeq) - {av^)z^2 {n^ “ ?T-eq) , (A.6) 
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with 


/ \ _ rjTself 

((T'u)self = 2 ^ 2 


re: 


(cJ'u)semi — 


Tsemi 


eq 


re, 


and {av^) 3^2 = 


73-J-2 


eq 


re, 


(A.7) 


eq 


Notice that one can also define 


which satisfies 


{(7v)2^3 = 


{av^)z-,2 = 
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Kinetic Equilibrium 

When self or semi-annihilations are non-negligible, because of crossing symmetry, the tran¬ 
sition amplitudes leading to them also induce scattering processes and these keep DM in 
kinetic equilibrium with its decaying products. In fact, when the DM particles become non- 
relativistic, the scattering processes occur at a faster rate than the annihilations because 


Dkin rejion-DM(re'u)kin re-non-DM 

- ^ - rvj - 

r ann reDM(o-u) 

ann rej^M 


(A.IO) 


where in the last step we take relativistic non-DM particles, which are more abundant during 
the radiation-dominated era. Because of this, the kinetic equilibrium is still maintained much 
after the DM freeze-out. 

In this work, when self and semi-annihilations can not be neglected, we further assume 
that all the non-DM particles are in equilibrium with each other. In particular, this means that 
the GB and the SM particles are in equilibrium. In order to avoid non-negligible contributions 
to Aeff) as explained in section 3.1, such equilibrium can only exist at temperatures greater 
than the muon mass. Because of this and because the freeze-out temperature is roughly 
mx/25, in this work we only consider self- and semi-annihilations in the GeV range. 

The situation is different in the opposite case, that is, when the 3—^2 mechanism 
dominates. In this case, non-DM particles do not participate in the DM production, and 
therefore scattering processes between DM and GBs or SM particles do not necessarily take 
place. Nevertheless, the 3 —)• 2 processes convert a significant fraction of the mass into DM 
kinetic energy. If the latter is not radiated away, DM heats up altering significantly structure 
formation [64-66] We can circumvent this by allowing sufficiently large DM couplings 
to the SM plasma or the GBs, so that the former remains in equilibrium with one of the 
latter [25]. As shown in section 5, the first possibility requires huge couplings because the 
scattering rate is suppressed by the electron Yukawa coupling. Therefore, we assume the only 
possibility left, that is, that GBs and DM are in equilibrium. Because of Eq. (A.IO), this is 
not in contradiction with the fact that self- or semi-annihilations are negligible. Notice also 
that in this case the DM mass is not required to be in the GeV range because we assume no 
equilibrium between the GBs and SM plasma. 

^ 2 (and 4 —>■ 2) processes are also constrained by the lifetime of the DM halo [67], however this 
constrain is typically subdominant. 
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B The Non-relativistic Approximation 


In order to calculate 7 for the different processes with only DM particles in the initial state, 
we take = exp(—ill/r) in Eq. (A.3). For temperatures much smaller than the DM mass, 
because of the Boltzmann suppression, the integration over the phase-space of the initial state 
is dominated by small three-momenta. In fact, by taking each of them equal to zero one finds 




^eq 

Amx 





d^Pj \ 

2Ej(27r)3 I 


(27r)‘^(5‘^(pa 


P/ 3 ) X 
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where the subscript NR means that all the initial state particles are taken at rest. After per¬ 
forming the two-body final state phase-space integration, this implies that the cross-sections 
are 
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(B.4) 


Notice that for the self- and semi-annihilation cases, these expressions are just the s-save piece 
of the partial wave expansion of the thermal averaged cross-sections. 

In Table 2, we show all the square amplitudes that contribute to these cross-sections 
when Xhx^ X^px 6 are neglected. For the sake of completeness, we also include the square 
amplitude for other processes with three DM particles in the initial state: XXX —?■ Xr] and 
XXX —>■ ijr]. Notice that when the 3 —)• 2 mechanism is relevant - that is, when either Z or 
Ajv is order one and A3 is small - those processes have a suppressed amplitude. We also show 
some processes with four DM particles in the initial state. Because Ax and Z are bounded 
from perturbative unitarity considerations as implied by Eq. (5.4), those processes can not 
dominate over the 3 —)• 2 mechanism. 
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